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Coherent exciton transport in 
semiconductors 



1.1 Introduction 

An exciton is a particle-like neutral excitation of solids and molecules composed of one 
electron and one hole bound together by the mutual electrical attraction [13411206112211 
122811286] . Its creation through internal charge separation is most frequently caused by 
the absorption of light and its demise is occasioned by electron-hole recombination, 
mostly with emission of light and less frequently non-radiatively. The many-electron 
ground state of the system, being an insulator, is immune to excitation until the 
excitation energy reaches a threshold G known as the energy gap. When external 
influences such as the electromagnetic field and lattice vibrations arc ignored, the 
exciton may be viewed as a robust state of an excited electron plus the hole which 
has been left behind in the valence electron states |230j . The hole acquires its positive 
charge from the loss of an electronic charge from the ground state whose total charge 
is neutralized by that of the ions in the molecule or solid. 

The photon-exciton interaction is responsible for the optical excitation (though 
not necessarily in the visible frequency range) of the exciton and for its spontaneous 
recombination emitting a photon (a quantum unit of light). The dipole matrix element 
responsible for the transition between the energy states is strong when the electron 
and hole wave functions overlap in space or match in wave vector. From Planck's law, 
the frequency of the emitting light Ex /h is proportional to the energy loss Ex in 
returning the exciton state back to the ground state, with h being Planck's constant. 

If the constituent electron and hole of the exciton are mostly localized at an ion, 
the exciton is localized, but with some probability to hop from site to site. Such a 
Frenkel exciton is common in molecules and molecular solids. At the other extreme, if 
the electron and hole wave functions are widespread as extended orbitals in a molecule 
or Bloch waves in a crystal, their bound state as the exciton can have their center of 
mass moving through the system with ease. Such Wannier (or Wannier-Mott) excitons 
are most common in broad-band and small-gap semiconductors (a semiconductor is 
distinguished from an insulator qualitatively by a smaller energy gap, with the fre- 
quency of the emitting light from the exciton spanning the range from visible light to 
very far infrared). Wannier excitons resemble the hydrogen atom or, more closely, the 
positronium system composed of an electron and a positron. Because of the dielectric 
screening of the electrical force in small-gap materials and sometimes the small effec- 
tive mass of the electron, the Wannier exciton radius is 10 to 100 times larger than 
the positronium radius, which is approximately 0.1 nm. 
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Excitons, being made of two fermions, behave as bosons on the scale larger than the 
exciton radius and therefore may macroscopically occupy a single quantum state [201 
I175ll40lll25ll3]|124ll77] . If the exciton lifetime is long enough to allow for reaching quasi- 
equilibrium, the dilute and cold gas of optically generated excitons may undergo Bose- 
Einstein condensation (BEC) [50l[87lfT98] . The critical temperature for exciton BEC, 
of the order of 1 K for typical densities in semiconductors, is basically the temperature 
at which the thermal De Broglie wavelength becomes comparable to the average inter- 
exciton separation. The possibility of achieving BEC of excitons by shining light on 
solids has been thoroughly investigated in the last fifty years (see the reviews j99[|178i 
[lailSainellllllESaiSOlETlE^ semiconductors are particularly appealing 

for this goal as they may provide excitons with a a lifetime (hundreds of ns in bilayer 
structures [30j ) longer than the time required for cooling. 

In an indirect-gap semiconductor such as silicon, where the momentum of the 
exciton docs not match that of the photon, the excitons are generally formed after 
relaxation of optical excitations with initial energy much higher than the gap. The 
indirect exciton has a long lifetime because its recombination with the emission of a 
photon requires the conservation of momentum to be satisfied by the assistance of a 
lattice vibration or trapping by a defect. Consequently, the excitons have time to form 
a large pool known as an electron- hole drop (see the reviews |209ill22| ). Alternatively, 
the delay in optical recombination may be due to the symmetry of the crystal, as in the 
direct-gap oxide CU2O which has conduction and valence bands of like parity hence 
the optical dipolar transition is forbidden |110| . The chapter by Kuwata-Gonokami in 
Volume 1 focuses on the aspects of Bose- Einstein condensation of optically generated 
excitons in semiconductors. 

In a direct-gap material, the spatial separation of the electron and hole can be en- 
forced by housing them in two layers sufficiently close to maintain their electric attrac- 
tion j232 pi60| [7T| . The recombination of such indirect excitons may then be controlled 
by changing the electron and hole wave function overlap with an electric field [51 1300] . 
An interesting phenomenon is the laser spot excitation of these indirect excitons, lead- 
ing to the formation of two concentric luminous circles centered at the laser spot plus 
other localized bright spots randomly placed between the circles. Whereas the forma- 
tion of the inner ring [34[ is due to the migration of indirect excitons away from the 
laser spot as optically inactive excitons, the localized bright spots [34] as well as the 
outer ring [Mll238l[571l204| form on the boundaries between electron-rich and hole-rich 
regions. At low temperature the outer ring is a necklace of evenly spaced bright spots, 
whose origin is not fully understood I103j . This system of excitons in a double quantum 
well is considered a good candidate for condensation. These experiments and related 
work in double quantum wells are reviewed in |2371[30ll31] l32 1l256llll2j . 

Another possible — and elusive — mechanism of condensation of excitons as bosons is 
that excitons form spontaneously at thermodynamic equilibrium even in the absence 
of an optical excitation. Such process signals the transition to a permanent phase 
known as excitonic insulator (EI) [T77 l [T34 l [r23 [ l55 [ [TT4] . which is originated by the 
instability of the normal ground state of either a semiconductor or a semimetal against 
the spontaneous formation of bound electron-hole pairs. The wave function of the 
strongly correlated EI ground state is formally similar to that proposed by Bardeen, 
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Cooper, and Schrieffer for superconductors [T^. As a matter of fact, both cxcitons 
and Cooper pairs are absent except as fluctuations in the normal high-temperature 
phase and form only in the ordered, low-temperature phase — respectively the EI and 
the superconductor. Besides, both condensation of excitons and that of Cooper pairs 
are best described in the reciprocal space of the crystal solid. The EI phase is reviewed 
in [Ml[n51iggi[5ni[T551[571 [Wl[T751[^ll56llll2[[M]] . 

It is intriguing to observe that condensation of other types of bosons composed of 
two fermions leads to spectacular manifestations of quantum mechanical coherence, 
such as the superfluidity ensuing from the pairing of "^He atoms [149] . Fermi alkali 
atoms confined in optical traps [211180], nucleons in neutron stars [79lll96| . the super- 
conductivity induced by Cooper pairs in metals [54] . and the non classical momenta 
of inertia in nuclei |1691 l23j. The above phenomena may regarded as distinct realiza- 
tions of superfluidity, associated to the coherent, dissipationless flow of charge and / 
or mass. However, excitons are neutral and stay dark unless recombine radiatively, as 
shown in Table [O] which compares the distinct features of the condensates made of 
composite bosons. The signature of the macroscopic order of the exciton condensate 
is, at present, controversial for the superfluid transport but its other manifestations 
will be discussed next. 

The aim of this Chapter is to illustrate some recent theoretical proposals con- 
cerning the detection of coherent exciton flow j2121l213ll214] . The reader may refer to 
the literature reviewed in Sec. 11.21 for a discussion of the conceptual and experimen- 
tal difficulties inherent in the realization of exciton condensates. Here we set aside 
such difficulties and adopt in a pedagogical way the simplest mean-field description 
of the condensate, on which we lay our theoretical development in order to detect the 
transport properties of the exciton condensate. 

In particular, we focus on the exciton analogues of two phenomena, i.e., Andreev 
reflection and Joscphson effect, which arc hallmarks of superconducting behavior, and 
stress the crucial differences between excitons and Cooper pairs. Our first main conclu- 
sion is that the excitonic insulator is the perfect insulator in terms of both charge and 
heat transport, with an unusually high resistance at the interface with a semimetal — 
the normal phase of the condensed state. Such behavior, which should be contrasted 
with the high electrical conductance of the junction between superconductor and nor- 
mal metal, may be explained in terms of the coherence induced into the semimetal by 
the proximity of the exciton condensate. Then we show that the exciton superflow may 
be directly probed in the case that excitons are optically pumped in a double-layer 
semiconductor heterostructure: we propose a correlated photon counting experiment 
for coupled electrostatic exciton traps which is a variation of Young's double-slit ex- 
periment. 

We last mention that, due to the interaction between electrons and light, not only 
can an exciton decay irreversibily into a photon or vice versa, but it can also exchange 
roles with the photon in a quantum-mechanically coherent fashion. Thus, the exciton 
may exist in the solid in the superposition state of an exciton and a photon, known as 
polariton. Whereas the photon energy varies linearly with its momentum at the speed 
of light in the vacuum, the exciton energy depends on the square of its center-of-mass 
momentum. For small momenta, the exciton and the photon can approximately match 
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Table 1.1 Excitonic insulator (EI) versus superconductor a la Bardeen-Cooper-Schrieffer (BCS). 
The interface referred to in the Table is the junction between normal and condensed phase. For a 
general discussion of the condensates made of composite bosons see [136j . For specific EI features 
see [114) (Meissner effect), |297j (superconductivity), [298j (superthermal conductivity), [21211213] 
(Andreev reflection), and [214| (Josephson oscillations). 



J. ilVOlLcli UiUJ-^v-'i by 


lllJ\.\^L LrUiiiU lilo LLld bKJL 


RiiS-lik'P QiiTiAT'r'nnrinpl'nT 


Nature of the composite boson 


Exciton 


Cooper pair 


Boson charge 


Neutral 


2e 


Boson momentum 


Crystal momentum 


Crystal momentum 






(commonly ignored in the free 






electron gas approximation) 


Boson mass 


Effective mass 


Effective mass of the electron 






quasiparticle in the Fermi level 






region (of thickness provided by 






phonon Debye frequency) 


Type of long-range order 


Diagonal 


Off-diagonal 


Superfluidity 


? 


Superconductivity 


Meissner effect 


No 


Yes 


Superthermal conductivity 


No 


No 


Nature of the quasiparticle 


Electron (hole) 


Bogoliubon 


Andreev reflection 


Yes 


Yes 


Interface electric conductance 


Decreased 


Increased 


Interface thermal conductance 


Decreased 


Decreased 


Proximity effect 


Yes 


Yes 


Josephson oscillations 


Yes 


Yes 



both their momentum and energy values, the coupling mixing the two states into two 
superpositions of photon and exciton with an energy splitting. Thus, the massless 
photon is slowed down by the massive exciton by virtue of the quantum-mechanical 
superposition. The Chapter by Yamamoto in Volume 1 deals with aspects of polariton 
condensation. 

The structure of this Chapter is the following: After a review of previous work 
(Sec. I1.2p . in Sec. 11.31 we illustrate the mean-field theory of the EI emphasizing its 
relation with the BCS theory of superconductors. We then introduce the phenomenon 
of Andreev reflection in Sec. 11.41 and analyze its observable consequences in Sec. 11.51 
Section [L6l on the Josephson effect ends the Chapter. 
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1.2 Physical systems 

This section briefly reviews recent theoretical and experimental works on exciton con- 
densation, focusing on diverse physical systems. Without attempting an exhaustive 
review, we refcr the reader to more comprehensive essays whenever available. 

1.2.1 Bose-Einstein condensation of optically generated excitons 

The pursuit of Bose-Einstein condensation of optically generated excitons in semi- 
conductors, which dates back to the sixties, presently focuses on both classic systems 
such as CU2O and novel low dimensional structures (for reviews see [99, 178,87, 1571 
[TM[I12[23ll[120l|30l[3ll[T5i|3a^ A very active field concerns 

"indirect" excitons. Such excitons are made of spatially separated electrons and holes, 
hosted in two quantum wells that are sufficiently close to maintain electrical attrac- 
tion between the carriers of opposite charge. This setup has several advantages: (i) 
The overlap of electron and hole wave functions is controlled by applying an elec- 
tric field along the growth direction of the bilayer heterostructure, thus increasing 
the exciton recombination time by orders of magnitude with respect to the single- 
well value [5l l300j . (ii) The confinement effect along the growth direction increases 
the exciton-phonon scattering rate, improving exciton thermalization [292^ . (iii) The 
dipolar repulsion among indirect excitons disfavors the formation of biexcitons and 
electron-hole droplets |38II15 6"246"222"267':148':269B7] as well as effectively screens the 
in-plane disorder potential |11 1. 220, 102, 107.208. 7]. (iv) As the electric field parallel to 
the growth direction may be laterally varied using suitably located electrodes, one may 
tailor the in-plane effective potentials for excitons, thus realizing artificially controlled 
traps P^l9 7,44,86, 75. 102, 107,223, 104, 7,8.. ramps [MfTi] . lattices [295,296,208,207^, 
"exciton circuits" [flJTl[TU51 IS5] . and "exciton conveyers" [175]. 

Exciton traps may also be created by means of the uncontrolled in-plane disor- 
der of the double quantum well j299| [36 l [34 lll02|ll07] . the strain experienced by the 
heterostructure [2571[Illll83l[179l[268l[2H5] , the laser-induced confinement [98l[9], the 
magnetic field |46| . The realization and control of exciton traps is a key capability to 
reach exciton BEC: As the long range order in two dimensions is smeared by quan- 
tum fluctuations, a weaker requirement for the macroscopic occupation of the lowest 
exciton level is that the exciton coherence length exceeds the trap size [31] . 

The present evidence of exciton BEC is based on distinct features of the emit- 
ted light (photoluminescence, PL) that appear at low temperature: (i) The PL dy- 
namics exhibits bosonic stimulation of the scattering of hot optically dark excitons 
into optically active low-energy states [3S]. (ii) The PL signal becomes noisy in a 
broad range of frequencies, as it occurs in the presence of coherence [38[|139[ri42| . 
(iii) The degree of polarization of the emitted light increases with decreasing tem- 
perature [144[ll45[ll03j . consistently with gauge symmetry breaking, (iv) The exciton 
mobility is enhanced, which may be attributed to superfluid behavior [33]. (v) The 
radiative decay rate increases, which may be explained in terms of "superradiance" of 
a macroscopic dipole [33] or collective behavior at the onset of condensation |144lll45j . 
(vi) The PL lineshape narrows and departs from the Maxwell-Boltzmann distribu- 
tion |139(ll44lll45lll46] . as it may be expected for the macroscopic population of a 
single exciton state. 
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However, some of the signatures |110|I239|I155| [72 ] listed above, taken separately, 
may have different explanations than exciton condensation |187|ll88llll8j (for a dis- 
cussion see Ref. [3T]). The most compelling evidence of BEC is probably the direct 
measure of coherence through interferometric techniques |2821I256||1031I1041I227II1051 I5] , 
which accesses the macroscopic exciton wave function in real space. On the theory 
side, the light emitted by excitons just after the onset of condensation is predicted to 
be coherent [1901 [Ml HH], with a sharply focused peak of radiation in the direction 
normal to the quantum- well plane |12HI294j . Besides, the instability leading to the 
external ring of evenly placed bright spots discussed in the introduction [311137] is 
possibly linked to exciton quantum degeneracy |153| . 

An intriguing issue is the role played by spin |165j in exciton condensation. In bilay- 
ers, the exciton spin is the component Jz of the angular momentum along the direction 
perpendicular to the planes, discriminating between optically active (Jz = zLh) and 
inactive states (J^ — zt2h). The most urgent questions concern the multi-component 
nature of the condensate [6711283] , the possibility of dark-exciton condensation [48ll49] , 
the role of spin-orbit coupling [95l l39l [281112331 . This research is fueled by the recent 
experimental evidence that the spin-relaxation time of indirect exciton is long and con- 
sequently exciton spin transport is long-ranged |152|I138] , as well as that spin textures 
and polarization vortices appear together with the onset of long-range coherence jl03| . 

1.2.2 Excitonic insulator in mixed- valence semiconductors 

In principle, any intrinsic semiconductor that may be turned into a semimetal, either 
by applying stress or by suitable alloying, may undergo a transition to the permanent 
EI phase. Favorable conditions are the presence of an indirect gap, which weakens the 
detrimental effect of dielectric screening on the exciton binding, as well as the nesting of 
electron and hole Fermi surfaces, which maximizes electron-hole pairing. Nevertheless, 
early experiments focusing on simple materials, such as divalent fee metals (Ca, Sr, 
Yb) and group V semimentals (As, Sb, Bi) were unable to confirm the existence of the 
EI. References [Ml[l35l[Ml[5Ql[ll6l[87l[l57l[l76l[l5^ review the work on the EI. 

Recently, a few experiments have pointed to the realization of the EI phase in 
mixed- valent semiconductors. The first class of candidate materials consists in rare- 
earth chalcogenides, such as TmSe^jTei-^; |28I272| . Smi-^jLa^jS j272l270j . Smi_a;Tma;S, 
YbO and YbS |270j . These compounds all crystallize in the NaCl structure and undergo 
a semiconductor-semimetal transition as the band gap G is changed from positive to 
negative values by applying high hydrostatic pressure to the sample. 

When the direct gap of TmSeo.45Teo.55, formed between the localized 4/^'^ levels 
and the 5d conduction-band states, is closing with external pressure, an indirect band 
gap develops between the highest valence Tm 4/^^ level ri5 at the F point and the 
mimimum of the A2' conduction band 5d states at the X point of the Brillouin zone. 
As the otherwise localized 4/ band is broadened and shows a maximum at F due to 
p(Se,Te)-/(Tm) covalent hybridization [113| . it is tempting to use a simple two-band 
model for interpretation, similar to the one ilustrated in Sec. 11.31 On the basis of low- 
temperature resistivity and Hall mobility measurements, the authors of [28] attribute 
the resistivity increase with the vanishing gap to a condensation of free carriers into 
excitons, placing the EI phase between semimetal and semiconductor, close to G ~ 0. 
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Later, the same group has reported a hnear increase of thermal conductivity and 
difFusivity with decreasing temperature and attributed it to exciton superfluidity [271] . 
Fehske and coworkers [26 ,27, 289] have suggested theoretically that the EI phase in the 
pressure-temperature phase diagram is narrower that the experimental claim, being 
surrounded by a "halo" regions made of preformed excitons coexisting with the normal 
semiconductor phase. The presence of this halo, precursor of the EI, explains the 
experimental findings and rules out the idea of a heat supercurrent, which conflicts 
with the general argument |1141I298| that a flowing condensate carries no entropy and 
thus no heat. 

Other candidate systems for the EI phase are the transition-metal chalcogenides 
TiSez [2T7l[mi[258l[THm[M7l [Tnl [TMl[ml l4^ 1 73l[T67l l42] . TazNiSes [274] . 

TaSe2 Eig, and the possibly ferromagnetic EI Gdl2 [THJ[153] (see Ref. HH] for a 
review) . The main evidence relies on the hole quasiparticle band structure, as extracted 
from angular-resolved photoemission [2581 [TBI [2471 [T0l [TMlfm] l4l [ml fTfllT^ . 
The much studied TiSe2, at a critical temperature of around 200 K, develops a charge 
density wave [217] which does not fit the standard model based on Fermi surface 
nesting [9011216] but it is consistent with the presence of an EI (see Sec. II. 3|) . In fact, 
the spanning wave vector of the charge density wave is the distance in reciprocal space 
between Ti 3d-electrons and Se 4p-holes, which are bound by Coulomb attraction. 
Therefore, the excitonic instability drives the charge density wave and may possibly 
couple with a periodic lattice distortion |247pi95p43pi74pi73pi67l[42] . though ahernate 
scenarios |13ip264ll216j have been suggested. Recent time-resolved photoemission data 
link the artificially induced collapse of the charge-ordered TiSe2 state to screening due 
to transient generation of free charge carriers, supporting the excitonic origin of the 
phase transition [21011100] . 

A third class of candidate systems consists in Kondo insulators [2] and heavy- 
fcrmion materials |147| . which arc mixed- valence semiconductors characterized by a 
flat /-type valence band plus a dispersive — say d-type — conduction band, typically 
exhibiting strongly correlated behavior. Such systems (e.g. SmBg) are often modeled 
by the Falicov-Kimball Hamiltonian, which takes into account the strong inter-band 
Coulomb interaction [66ll200p 21 H [2T5]. Sham and coworkers have shown [20T][200] that 
the exciton condensate made of / holes and d electrons may spontaneously break the 
lattice inversion symmetry and lead to a ferroelectric phase transition of electronic 
origin, whereas conventional ferroelectricity is associated to lattice distortion [133] . 
The predicted experimental signatures, supported by some evidence [273[[83] . include 
the divergence of the static dielectric constant, a ferroelectric resonance in the mi- 
crowave absorption spectrum, and a nonvanishing susceptibility for second-harmonic 
generation. 

If intraband hybridization dominates over Coulomb interaction, then the exciton 
condensate wave function acquires a different type of symmetry — ^p-wave — which ex- 
cludes the ferroelectric scenario but allows the coupling with the lattice [60]. In this 
latter case the excitonic instability manifests itself as a spontaneous lattice deformation 
which may explain some of the phase transitions known as ferroelastic [29] . 
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1.2.3 Permanent exciton condensation in bilayers 

In order to investigate permanent exciton eondensation in semiconductor bilayers, one 
strategy is to host electrons in the first layer and holes in the second layer [23211160) . 
This task is nowadays accomplished by means of suitable electric gates which allow 
to separately contact the layers j235lfTT6l[T99lfT27ll2021l251j . The spacer beetween the 
two quantum wells suppresses the inter-layer tunneling which induces exciton recom- 
bination, but it is sufficiently thin to provide strong inter-layer Coulomb interaction 
(see Ref. [92] for a recent review). This setup allows for measuring the Coulomb drag 
resistance, which is the inverse ratio of the electric current measured in one layer to 
the open-circuit voltage developed in the other layer in turn. Such drag resistance is 
predicted to diverge in the presence of exciton condensation, as the exciton binding 
correlates the motion of carriers in the two layers .266. 1108llll5| . Recent measure- 
ments [52 ll2251 [5T| point to low-temperature anomalies in the Coulomb drag which 
may originate from an excitonic instability, though other strongly correlated phases 
are possible [92] . 

An alternate strategy is to place electrons in both layers in the presence of the 
magnetic field (see [63 ] 164 ] [259] for reviews). The field bends classical electron trajec- 
tories into circular cyclotron orbits. As such orbits may be placed all across the plane, 
overall their quantized energies consist in highly degenerate "Landau levels". Since 
the level degeneracy is the number of quanta of magnetic fiux that cross the plane, 
for sufficiently high fields and identical layers the lowest Landau level in each layer 
will be half filled by electrons (single layer filling factor ly = 1/2, total filling factor 
i^T — !)• Note that in this quantum Hall effect regime, routinely detected through the 
quantization of the Hall resistance. Landau levels may be considered cither half filled 
or half empty. Therefore, one may switch to the excitonic parlance [14011284] . regarding 
one layer as filled by electrons and the other one by holes: In this picture the exciton 
"vacuum" has the lowest Landau level totally filled in one layer (v = 1) and empty in 
the other layer (i' = 0), thus excitons are created by moving electrons from one layer, 
which leaves a hole behind, to the other one [164) . 

There is significant evidence, based on low-temperature transport experiments, 
that the bilayer ground state is a condensate of excitons. The first hint is a huge 
enhancement of inter-layer tunneling solely due to many-body effects, clearly pointing 
to strong inter-layer coherence [244) . The most compelling observations are based on 
counterfiow measurements [232)1160)1248) , where the electric currents of opposite sign 
and like magnitude that flow in the two layers provide zero total electric current and 
a net exciton flow. For fllling factors other than = 1 the Hall voltages separately 
measured in the two layers are equal and opposite in sign, whereas for z^t = 1 they 
both drop to zero, consistently with the flow of an uncharged object such an exciton 
[126)1260)1277) . To prevent edge states — always present at the boundary of quantum 
Hall systems and unrelated to excitons — from playing a role in transport, the Coulomb 
drag has been recently measured in the "Corbino" annular geometry, confirming the 
excitonic nature of transport [2551l254[[69lll80) . whereas the superfiuid character of 
the exciton flow is unclear. The above scenario is supported by the measurement of 
quasiparticle and collective excitations by means of tunneling [245) [8T] and inelastic 
light scattering [161)1117] spectroscopies. 
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Interesting theoretical predictions concern the response of the bilayer exciton con- 
densate to external electromagnetic fields [1^1115 [[2151152] and impurities [HI], as well 
as the transport properties of hybrid circuits including exciton condensates and super- 
conductors [57 lll92ir242j . For weak inter-layer interaction or filling factors other than 
i^T = 1 , bilayers are predicted to undergo phase transitions to other strongly correlated 
phases, such as paired two-dimensional Laughlin liquids and Wigner solids [284], or 
peculiar excitonic charge density waves |45j . 

1.2.4 Graphene-based systems 

Graphene — a recently discovered allotrope of sp^ bonded carbon — is a one-atom thick 
two-dimensional honeycomb lattice [76 lll85ir219[ll37j . Its peculiar electrical and me- 
chanical properties — chemical stability, high mobility, easiness of making electric contacts — 
have stimulated observations by means of different electron spectroscopics and scan- 
ning probes. Intensive investigations have uncovered new physics (e.g. Klein tunnelling, 
anomalous types of quantum Hall effect), rooted in the unusual character of quasipar- 
ticle excitations, that, in the neighborhood of the Fermi energy, are massless chiral 
Dirac fermions. In fact, conduction and valence bands form specular cones whose 
apexes touch in the two inequivalent points K and K', located at the corners of the 
hexagonal two-dimensional Brillouin zone. These two points, which map into each 
other by a rotation of 27r/6 [51], are the Fermi surface of the undoped system, hence 
graphene is a zero-overlap semimetal. 

In principle, graphene is a good candidate system for EI, since: (i) the density of 
states vanishes at the charge neutrality point, hence the long-range Coulomb interac- 
tion is unscreened (ii) the perfect electron-hole symmetry of Dirac cones favors the 
nesting of electron and hole isoenergetic surfaces. Khveshchenko |130) was the first to 
suggest that graphene hides a latent excitonic insulator instability. The EI phase is a 
charge density wave alternating between the two inequivalent triangular sublattices, 
its spamming wave vector connecting K and K' in reciprocal space. A stack of graphite 
layers in a staggered (ABAB...) configuration, with the atoms located in the centers 
and corners of the hexagons in two adjacent layers, respectively, could stabilize the EI 
by enforcing interlayer Coulomb interaction. 

After this seminal prediction, many theoretical works have tried to estimate the size 
of the EI transport gap as well as the stability of the EI phase (see for example |262[ 
[58]r73]l84]l275] and references therein as well as the reviews [I84lll37p226j ). The absence 
of consensus is not surprising, as the many-body problem in graphene is presently an 
open issue |137[|261j . Experiments show that electrons in graphene allegedly behave 
as non-interacting particles [76lll85i[219) . except for small effects related to velocity 
renormalization [65] , coupling with phonons [76 lll85|[219j / plasmons [25l (here we are 
not concerned with the fractional quantum Hall effect [59ll24] , induced by the magnetic 
field). Therefore, if the EI energy gap ever exists, it must be smaller than the present 
spectroscopic resolution. 

A related theoretical proposal concerns permanent exciton condensation in double- 
layer graphene |159[ [56 ] . The idea is to separately contact the two layers, which are 
spaced by a dielectric medium, in order to induce the same quantity of charge with 
opposite sign in the two layers. With respect to the double-layer made of usual semi- 
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conductors mentioned in Sec. I1.2.3| here the advantage is the smaller value of the 
transverse electric field required to polarize the bilayer, due to the zero energy gap of 
graphcne. The estimate of the Kosterlitz-Thouless temperature required to undergo 
the EI phase is debated theoretically [T7nl[Mn[ml [Tm[78l [T2l[T7a[T58l[MTl[25^ : 
recent Coulomb drag measurements |132l[85] point to the importance of inter-layer 
interactions. 

In the absence of a dielectric spacer, undoped bilayer graphene is predicted — among 
other proposals — to undergo an excitonic ferroelectric phase that spontaneously breaks 
which- layer symmetry and polarizes the layers in charge |171[|181ll291| . The excitonic 
instability, which opens an energy gap, appears to sensitively depend on the interaction 
range [9Tj: in the case of finite range, the expected electronic phase is nematic and 
gapless |151[|263) (see reviews jl631ll82j ). The experimental observation of a transport 
gap at the charge neutrality point is controversial |276[|168ll265l[7D] . 

1.3 Two-band versus BCS model 

In this section we contrast the mean-field theory of the EI to the BCS theory of 
superconductors. Wc compare the equations of two model junctions, (i) one between EI 
and scmimctal (SM-EI), (ii) the other one between superconductor and normal metal 
(N-S). In both cases the phase boundary is due to the variation of the order parameter 
that changes along the direction perpendicular to the interface, tending respectively to 
zero in the bulk normal phase and to a constant value inside the condensed phase. The 
quasiparticle amplitudes for both SM-EI and N-S junctions are formally identical and 
are used in Section [1.41 to compute the flow of charge and heat through the junction. 

1.3.1 The SM-EI junction 

We start studying the junction between semimetal (SM) and EI on the basis of a 
spinless two-band model. The SM has overlapping isotropic conduction and valence 
bands {b and a, respectively) of opposite curvature and one electron per unit cell, 
hence the Fermi surface is a sphere in momentum space, located at the nesting of the 
two bands; since there are TV electrons and 2N states available, the nesting occurs at 
zero energy in Fig. 11.11 One may turn the SM into an EI by either changing the SM 
stoichiometric composition through suitable alloying or applying stress, which opens 
a gap of size 2A in the bulk EI in virtue of the strong intcr-band Coulomb interaction 
[cf. right panel of Fig. Il.lf a)]. The variation of the EI order parameter A(r) (defined 
below) along the coordinate z normal to the interface determines the effective interface 
potential, as shown in Fig. I1.2f a) . 

The Hamiltonian of the SM-EI junction is 



Here Hi is the kinetic term which embodies the effect of the ideal and frozen crystal 
lattice on electrons with the envelope function in the effective mass approximation 



'HsM-EI — Ho + Hi + ■^2- 



(1.1) 




i—a,b 



(1.2) 
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Fig. 1.1 Junction between semimetal (SM, left) and excitonic insulator (EI, right), (a) Quasi- 
particle energy to vs wave vector . The labels mark the allowed elastic scattering channels for 
an incoming electron (labeled /) with fcz = q'^ ■ In particular, A is the interband (Andreev) 
reflection, B the intraband reflection, C the intraband transmission, and D the interband 
transmission. The size of the EI gap is 2 A. (b) Isoenergetic contour lines in the {kx, kz) space 
for the energy uj shown in panel a. The arrows point to the group velocities of electrons in 
the different scattering channels. 

The field operator ipa{r) [ipb{r)] annihilates an eleetron in the valence (conduction) 
energy band at the position r in space. The real-space band operators Si{r) appearing 
in Eq. (|1.2p take the form: 

£a(r) = G/2 + (2m)-iV2; (1.3a) 
eb{r) = -G/2 - (2m)-^V^ (1.3b) 

Here m is the (positive) eiTectivc mass, G is the (positive) band overlap, and the 
energies are measured from the Fermi surface [114] . Throughout this work we put 
h = 1 and assume that the system has unit volume. The valence- and conduction- 
band energy levels of the non- interacting bulk crystal, eigenvalues of Jii, are 



ea{ka)=G/2-{2m)-'kl 
etikb) = -G/2 + {2m)-^kl 



(1.4a) 
(1.4b) 
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where ka and fef, refer to the respective band extrema. We assume the valence-band 
has a single maximum at fe = whereas the conduction-band a single minimum at 
k ~ w, and ignore complications due to the presence of equivalent extrema. The Fermi 
wave vector is given by fcp = mG. The two-body term 7^2 consists of the inter-band 
Coulomb interaction, 

% = jdrdr'i^l{r)i:l{r')V2{r-r')i',,{r')^a{r), (1.5) 

with V2(r) being the dielectrically screened Coulomb potential |123| . Renormalization 
effects due to intra-band Coulomb interaction and temperature dependence are taken 
into account into the energy band structure (|1.4p . The one-body term "Hg is the sum 
of two parts, 

Ho = V + Vhyb. (1.6) 

V is the intra-band term, 

V= J2 J<ir4ir)Vir)Mr), (1.7) 

i—a,b 

which includes the effects of the band offset as well as those of possible impurities and 
defects at the interface, such as a thin insulating layer, via the singlc-particlc potential 
V{r). The potential V{r) can also include the effect of a voltage bias applied to the 
junction in a steady-state regime. The inter-band term, 

Vhyb = Jdr^l{r)Vi,yi,ir)^Pair) + H.c, (1.8) 

hybridizes b and a bands by means of the potential V^iyb(i'). This term may be orig- 
inated e.g. by the change of an element in the SM compound. The influence of the 
potential Vhyb('') on exciton condensation, which depends on the symmetry of the 
bands involved, is discussed in Refs. |20HI200| l60]. 

Both the hybridization potential and the EI order parameter contribute additively 
to the small band gap formed between the two overlapping bands [cf. right panel in 
Fig. ll.lf a)] hence their separate contributions cannot be distinguished by spectroscopy, 
as it was early recognized |1141I136[|93] . However, in this Chapter wc show that the 
transport properties across the SM-EI junction are distinctive features of the EI. The 
key point is in the length scale of the variation of the effective interface potential 
which reflects or transmits the electron. A component of the effective potential is the 
position dependent order parameter A(r), which decreases from the bulk value in the 
EI region to zero in the SM region [see Fig. I1.2r a)]. The length scale of the change 
is the coherence length in the EI, much longer than the lattice constant. Since this 
scenario is most likely if the lattices of the two components are as similar as possible, 
we classify as homogeneous the junction with A(r) being the only contribution to the 
interface potential. 

On the other hand, the one-electron interface potential 14yb('") due to the change 
in hybridization has an abrupt variation on a length scale of the order of a few atomic 
layers, thus considered in our context as a heterogeneous junction (this includes the 
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case of a Schottky barrier). We have studied the common physical features of the 
heterogeneous junction, including the abrupt band edge discontinuity, the short-ranged 
interface potential, and the impurities at the interface. Whereas the charge carriers in 
the heterogeneous junction experience uninteresting intraband reflection, electrons in 
the homogeneous junction change valley when backscattered as a feature of the EI band 
mixing. In the rest of the Chapter we focus on the homogeneous junction and refer 
the reader elsewhere [212) for the study of the heterogeneous junction (corresponding 
to the case Vhyb('') 0)- 

Following Sham and Rice |230| . we introduce the electron quasi-particle amplitudes, 



f{r,t) 



{^o\Mr,t)\n) 



(1.9a) 
(1.9b) 



Here |4'o) and arc the exact interacting ground states with N and A^+1 electrons, 
respectively; the quantum index k labeling the electron quasiparticle means the crystal 
momentum only in the bulk phase as the overall translational symmetry is destroyed 
by the presence of the junction. States and operators are written in the Heisenberg 
representation [1] (flagged by the tilde symbol on operators): 



V;,(r, t) = cxp(i [n - fiU] t) il>,{r) exp(-i [U - t) , 



(1.10) 



where ^ is the chemical potential (here ^ 
number operator M is defined by 



due to electron-hole symmetry) and the 



E 



(1.11) 



Writing down the Heisenberg equations of motion for the operators tpi{r,t), exploit- 
ing the mean-field approximation to express them in a closed form, and neglecting 
unessential intra-band Hartree terms, we derive a set of two coupled integro-diffcrential 
equations for the amplitudes f{r,t) and g{r,t): 



.dfirvt) 

. dg{r,t) 
' dt 



_^_^+y(r) 

2m 2m 

Yl + ^ + Vir)\gi 
2m 2m 



f{r, t) + y dr'A(r, r') g{r', t) , (1.12a) 
r,t) + Jdr'A*{r',r) f{r',t) . (1.12b) 



The built-in coherence of the exciton condensate, A(r, r'), appearing in Eqs. (|1.12p 
for fc > fcp is defined as 



A(r,r') = V2{r - r') {^o\Mr) i^Hr') |*o) 



(1.13) 



Except for the factor V2{r — r'), A(r,r') is the exciton macroscopic wave function. 
In fact, (4'o|V'b('') i^Hf') l^o) is the average on the many-electron ground-state of the 
operator destroying an electron- hole pair, i.e., one b-band electron at r and one a-band 
hole at r' [the electron creation operator ipl^ir') may be regarded as a hole destruction 
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operator]. Such average is zero in the SM phase, since for k > k-p 6-band levels are 
empty and a-band levels filled, but it acquires a finite value in the EI phase. Besides, 
the finiteness of A(r,r') reflects the new periodicity in real space of the EI phase, as 
the electron density shows an additional charge-density-wave-like order characterized 
by the wave vector w displacing the extrema of a and b bands |135| . For k < kp the 
roles of electrons and holes are exchanged hence the definition ()1.13p of A(r,r') is 
modified accordingly. 

The built-in coherence A(r, r') generically depends on both ccnter-of-mass and 
relative-motion coordinates, but inside the EI bulk the center-of-mass part of the 
condensate wave function is a plane wave with zero momentum, hence A depends 
only on the relative coordinate i — r'. We expect A(r — r') to smoothly vanish when 
\r — r'\ becomes larger than the characteristic length, the exciton radius. This allows 
us to easily find the bulk solution of the system of eqns ()1.12alll.l2b[) [V{r) = 0] in 
terms of the two-component plane wave 

\9kir,t)J WJ 

with energy 



ioik) = ^ek + \^k\, (1.15) 

where = (fc^ — kp) /{2m) and A^ is the Fourier component of A(r). The amplitudes 
are normalized as 

^Ul + ^] , \ukf + \vkf^l, (1.16) 



(1.17) 



2 V Ek 

and the relative phase between Uk and Vk is given by 

Uk Afc 



Vk w(fc) - S.k 



When Afc = 0, the amplitude (|1.14l) is the trivial solution with Uk = I and Vk = 0, 
i.e., a conduction-band plane wave. When excitons form a condensate, solution (jl.l4p 
is admissible only if the self-consistency condition derived by the definition of A(r) 
is satisfied. This condition, which can be easily obtained from Eq. (|1.13p . is formally 
analogous to the BCS gap equation: 

with V2,fc being the Fourier component of V2{r). 

In general, the amplitudes f{r,t) and g{r,t) are the position space representation 
of the stationary electron-like elementary excitation across the whole junction. Taken 
together, they signify the wave function of the quasiparticle: / (g) is the probability 
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amplitude for an electron of belonging to the conduction (valence) band. They satisfy 
the normalization condition 



dr 



|/(r,i)l' + |5(r,0l'J =1, (1-19) 

and have always positive excitation energy ut due to the definitions (|1.9m.l0|) . The 
probability current density J{r,t) can be found starting from the definition of the 
probability density p{r,t) for finding either a conduction- or a valence-band electron 
at a particular time and place, p{r,t) = |/|^ -I- l^fj^. After some manipulation of the 
equations of motion (|1.12l) , one derives the continuity equation 

V-J = 0, (1.20) 



dp 
dt 



where 



J = lmlr-f-g*-g}. (1.21) 



m rn 



Note that the two terms appearing in the rhs of Eq. (|1.2ip . referring respectively to 
conduction and valence band electrons, have opposite sign since the curvature of the 
two bands is opposite. One can verify that the semiclassical group velocity of the 
quasiparticle, Vg = VfeW, coincides with the velocity v given by the full quantum 
mechanical expression (|1.2ip . with J — pv. 

1.3.2 The N-S junction 

In this section we highlight the suggestive parallelism between the formalism intro- 
duced in section (|1.3.ip and the treatment of quasiparticle excitations in conventional 
superconductors, as modeled by the BCS theory. 
The Hamiltonian of the N-S junction is 

nN-s = n'„ + n[+n'2. (1.22) 

Electrons in the metal experience the crystal lattice potential through H'^ , 

cr=t,4- 

The space-dependent field operator ipaif) annihilates an electron with spin a in the 
two-fold degenerate conduction energy band, whose energy dispersion is taken to be 
parabolic for simplicity, 

£{k) = k^/{2m), (1.24) 

e(fc) being the eigenvalue of H'l. With respect to the SM, the role of a and b bands 
is replaced by the two spin flavors f and J,. The Fermi wave vector fcp is fixed by the 
condition that there are N electrons in the system, with p — kp/2m. 
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Fig. 1.2 Andreev reflection, (a) SM-EI junction. An incoming fe-band electron is backscat- 
tered into tiie a band. The tliicker (tliinner) curves are tlie renormalized (bare) bands at 
different vafues of tlie z coordinate, witli A{z) being tlie corresponding order parameter. 
Only the relevant low-energy portion of the spectrum is shown here, with numbers from 1 to 
6 pointing to the time sequence of the reflection process (the arrows represent group veloci- 
ties). The inset illustrates that the reflection process may alternatively be seen as a coherent 
flow of excitons from the SM into the EI. (b) N-S junction. In contrast to panel a, an electron 
is backscattered as a hole, hence the whole process may be thought of as the dissipationless 
flow of Cooper pairs through the interface. 

The relevant two-body interaction is attractive and short-ranged, 

H'2 = -g ydr^j(r)7A|(r)Vi(r)Vt(r), (1.25) 

with g being a positive constant parametrizing the combined effect of Coulomb and 
electron-phonon interaction in the vicinity of the Fermi surface |54| . The short-range 
interaction (|1.25p does not affect electrons with parallel spin as a consequence of Pauli 
exclusion principle. The effective potential ()1.25p results from the competition between 
Coulomb repulsion and the screening effect of the positive ions in the lattice. Close to 



Two-band versus BCS model 17 



the resonance frequency of the ion motion, the ions give a very large response to the 
perturbation induced by an electron charge. The resulting cloud of the moving electron 
plus the polarized ions has a net positive charge, then inducing a weak electron-electron 
attraction, whose characteristic energy is a tiny fraction of the Fermi energy. As in 
the case of the SM-EI junction, the boundary between N and S phases is determined 
by the variation along z of the pair potential associated to 7^2, defined below [see 
Fig. I1.2f b)]. A residual effect of Coulomb interaction is to shift the energy levels, that 
are already renormalized in the dispersion relation (|1.24[) . The one-body term 



E 



dr^PUr)V'ir)Mr), 



(1.26) 



arises from the possible impurities and defects at the interface, as well as the applied 
bias voltage. 

In order to find out the quasiparticlcs of the N-S junction, we follow the same 
approach as for the SM-EI junction, with one important difference [12] that derives 
from the following definition of the amplitudes: 



9'{r,t) 



{%\i>l{r,t)\^i). 



(1.27a) 
(1.27b) 



Here IvP^) is the state with one quasiparticle added to the many-electron ground state 
Y^'q). According to eqns (|1.27p . the number of particles is not a constant of motion, 
as we add both an electron [eqn (|1.27ap ] and a hole [eqn (|1.27b[) ] to |^o). This is 
allowed within the gran canonical framework, where the chemical potential /i is the 
independent thermodynamic variable instead of N . 

The resulting 'Bogoliubov-de Gennes' equations of motion are: 



' dt 

. dg'{r,t) 
' dt 



2m 
2m 2m 



V'ir) 



/'(r,t) + A'(r)g'(r,i), 
g'ir,t) + A'*{r)f'{r,t), 



with the local pair potential A'(r) being defined as 

A'ir)^-g{%\^P^ir)rP^{r)\%] 



(1.28a) 
(1.28b) 

(1.29) 



The space-dependent parameter A'(r) may be regarded as the center-of-mass wave 
function of the condensate made of Cooper pairs. The latter are bound pairs of two 
electrons with opposite spins, as it is evident from the definition (|1.29p . The product 
of the two operators that destroys a Cooper pair, appearing in the so-called anomalous 
average (|1.29p . does not commute with the number operator as it breaks the U{1) 
gauge symmetry of total Hamiltonian "Hn-S- Nevertheless, since the number of electrons 
is macroscopic, the number fluctuations, of order VN, are small and may be neglected 
with respect to the average value of N. Moreover, since a value of one or two is still 
small with respect to Vn, the ground states with either N or N — 2 electrons must be 
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regarded as identical, so the anomalous character (i.e., breaking the gauge symmetry 
of "Hn-s) of the definitions (|1.29p and ()1.27p is physically irrelevant. 

It is remarkable to observe that the systems of equations (jl.l2p and (|1.28p . re- 
spectively describing the SM-EI and the N-S junctions, are formally identical in the 
homogeneous case. This corresponds to put respectively V{r) = in cqns (|1.12p and 
V'{r) = in cqns ()1.28|) . as well as to take the built-in exciton coherence in ()1.12|) as a 
local operator, A(r, r') = 6{r — r') A(r) (then Afc does not depend on fc). However, 
the quasiparticle amplitudes for the two model junctions signify profoundly different 
types of single-particle excitations. 

In the EI, to obtain a free electron in the 6-band — for k > — one has to break an 
exciton among those forming the condensate, that is a bound pair of a 6-band electron 
and a-band hole. The way to do this is to either create an electron in the b band, 
whose amplitude component is /, or destroy a hole in the a band, whose amplitude 
component is g. 

In the S, to obtain an unbounded single particle with spin f; one has to break 
a Cooper pair in the condensate. This is accomplished by means of either creating 
an electron with spin f or destroying an electron with spin J,. The latter option is 
equivalent to creating a hole with spin f, as a consequence of time-reversal symmetry. 
The components /' and g' are the amplitudes for the propagation of the electron and 
the hole, respectively. Table [01 compares the key features of the EI with those of the 
S, with regards to both the ground state and the quasiparticle excitations. 

In the following we are interested in comparing the SM-EI and N-S systems. In order 
to stress their formal analogy, hereafter we drop the prime symbol to label quantities 
referring to the N-S junction (as /, g, A, etc) and use the same notation in both cases. 
With this convention, formulae (|1.14| - [TT^T|) obtained for the SM-EI junction hold for 
the N-S junction, too. 

1.4 Andreev reflection at the interface between excitonic insulator 
and semimetal 

In this section we introduce the phenomenon of Andreev reflection at the N-S bound- 
ary as a paradigm to discuss the transport through the SM-EI junction. There are 
three qualitatively important results that are common to both systems: (i) all three 
cartesian components of the velocity change sign when the quasiparticle is reflected 
(Sec. ll.4!T|) (ii) the ratio of incident quasiparticles C(w) which are transmitted through 
the interface depends on the coherence factors of the condensate, being strongly sup- 
pressed close to the gap (Sec. I1.4.2p (iii) the condensate induces pairing on the normal 
side of the junction (proximity effect, Sec. I1.4.3|) . 

1.4.1 Velocity inversion at the interface 

The electrical transport across the N-S junction exhibits high conductance behavior 
at vanishing applied voltage bias. This evidence seems to contradict the fact that 
quasiparticle excitations are gapped in the S (see Fig. II. 3|) : quasiparticles in the bulk 
N approaching the junction with energy smaller than the gap, < a; < A, cannot 
penetrate into the bulk S. This paradox is solved by Andreev reflection, which is 
illustrated below. 
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Fig. 1.3 Junction between normal metal (N, left) and superconductor (S, right). The plot 
shows the quasiparticle energy uj vs wave vector kz in the two bulk phases. The labels mark 
the allowed elastic scattering channels for an incoming particle (labeled 7) with kz — ■ A 
is the Andreev reflection, B the normal reflection, C the normal transmission, and D the 
cross-branch transmission. Note that the energy of the particle is positive: on the N side the 
hole branch for \kz\ < fcp is obtained by inverting the energy of the portion of band filled 
with electrons (showed as a dashed curve) with respect to the Fermi surface. 

Consider A(z) to be a smooth complex increasing function of z, tending respectively 
to the asymptotic values zero when z — > — oo. inside the bulk N, and Aq when z — > +oo, 
inside the bulk S [Fig. ll.2r b)]. Following Andreev [T2], we note that the medium under 
consideration is completely homogeneous with an accuracy 2m |Ao| /fcp — a very small 
quantity in typical superconductors. Therefore, we seek a solution of eqns (|1.28p in 
the form 

/(r) = e'^-"-,7(r) , g{r) = e'^— x(r-) , (1.30) 

where n is some unit vector and ?/(r) and xi''') ^^re functions that vary slowly compared 
to e'''f Substituting (|1.30p in (|1.28p and neglecting higher derivatives of 77 and Xi 
we obtain 




(iupn • V + w) ri[r) - A(z) x(r) = 0, (T31a) 
(iwpn • V - cj) x(r) + A*(z) ?7(r) = 0, (1.31b) 

where v-p — kp/m. It is easy to find for z — !• ±00 the asymptotic form of the solutions 
of eqns (jl.3ip describing the reflection of the quasiparticle falling on the junction. 
When z — )■ — c» we put A{z) = 0. Then 




(1.32) 
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Fig. 1.4 Sketch of Andreev reflection for quasiparticles approaching the junction from the 
normal-phase side, (a) SM-EI junction. The reflected left-going electron is equivalent to a 
right-going hole, (b) N-S junction. The reflected left-going hole is equivalent to a right-going 
electron. 

where n • fei = co/vp, n • k2 = —ui/v-p] Ci and C2 are arbitrary constants. The first 
term on the rhs of eqn (|1.32p corresponds to an electron whose velocity v (or J) lies 
along n, and the second term to a hole whose velocity lies in the opposite direction to 
n (in fact uj/v-p <C k-p since 2m |Ao| ^ 1). If > 0, then the wave function (|1.32p 
describes an electron incident on the boundary and reflected as a hole on the N side; 
if Uz < 0, it describes an incident hole reflected as an electron. When z — )• +00 we put 
A(z) = Aq in eqn (|1.3ip . The solution describing the transmitted wave (Jz > 0) has 
for uj > IAqI the form 



,i>^/2 



where C3 is a constant, ip is the phase of Ag, 



lAol" for >0, 



n • = —v-p ^\/ Lu^ — |Ao| for < 0. 



(1.33) 

(1.34a) 
(1.34b) 



As expected, for uj < \Aq\ the functions 77 and x decay exponentially as 2 — > +00, 
hence the quasiparticle is prevented from entering the bulk S. However, all three carte- 
sian components of the velocity of the reflected particle change sign [see Fig. I1.4f b)]. 
This remarkable behavior, which is not due to interface roughness since we take the 
interface to be completely flat, is the key to explain the electric transport through 
the N-S junction. An electron with velocity v is Andreev-reflected into a hole with 
velocity —v which carries exactly the same current as the incident electron. In fact, in 
virtue of time-reversal invariance, the hole moving with velocity —v may be regarded 
as an electron moving with velocity v [Fig. ll.4f b)]. Therefore, we may understand the 
process of Andreev reflection as an electron above the Fermi surface forming a Cooper 
pair with another electron below the Fermi surface on the N side: such pair moves 
to the S side merging into the condensate, whereas the second electron leaves a hole 
behind in the N Fermi sea |287) . 
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The results obtained in this section hold also for the SM-EI junction, provided one 
links 77 and x components to the probability amplitudes of an electron of being in 
either 6 or a band, respectively, as illustrated in Fig. I1.4f a). Apparently, the reflection 
process at the SM-EI junction seems the usual reflection of an electron from the gap 
barrier. However, the complete reversal of the velocity vector suggests that the reflected 
electron may be regarded as an incoming hole with the same velocity as the incoming 
electron. The idea is that the overall reflection process may be thought of as the flow 
of electron-hole pairs — excitons — from the SM to the EI side, where they merge into 
the exciton condensate. Below we substantiate this alternate interpretation. 

1.4.2 Coherence factors in the transmission coefficients 

To proceed we specify the functional form of the interface potential, assuming that 
the excitonic coherence is a step function at the SM-EI interface, A(r) = A9{z), with 
A > 0. Moreover, we introduce a simple-minded model for the effect of disorder at the 
origin (e.g. an insulating layer) through the (5-potential V{r) = H 6{z). In the following 
we abandon the slowly varying amplitude approximation and look for solutions of the 
full eqns (|1.28p . requiring a larger number of scattering channels than those used in 

Sec, on 

Carriers coming from the bulk SM with energies slightly outside the EI gap have, 
say for the incident electron at /, two reflection channels, A and B, and two trans- 
mission channels at C and D (see Fig. II. ip . If the energy lies within the gap, only 
the two reflection channels are possible. Whereas the interface — by breaking the lat- 
tice translational symmetry — can in principle connect different parts of the Brillouin 
zone |229j , here the relevant regions of the wave vector space are the two valleys near 
the gaps in the bulk EI for those states with the same component of the wave vector 
parallel to the interface. 

We consider the elastic scattering at equilibrium, matching wave functions of the 
incident (/), transmitted (C and D), and reflected {A and B) states at the boundary, 
following the approach of Ref. [22] . In the bulk EI, there are a pair of magnitudes of 
k associated with uj, namely 

fc± = V2^\/kl/2m ± (w2 - A2)i/2. (1.35) 

The total degeneracy of relevant states for each uj is fourfold: ±k^. The two states 
±k~^ have a dominant conduction-band character, whereas the two states ±k~ are 
mainly valence-band states. Using the notation 

the wave functions degenerate in uj are 

*±..= Mo±-^^ ^ f"°)e±--, (1.37) 

\voJ \uoJ 

with the amplitudes uq , vq defined as 
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[U (a;2 - A2)i/2 

possibly extended in the complex manifold. With regards to the SM bulk, A = and 
the two possible magnitudes of the momentum q reduce to = [2'm{kp/2m ± u})]^^'^ , 
with wave functions 

vI/±,+ = Qe±'«^^ = Qe±'^"^ (1.39) 

for conduction and valence bands, respectively. The appropriate boundary conditions 
are: (i) Continuity of * at 2 = 0, so *ei(0) = «'sm(0) = ^(0). (ii) [/^i(0) - /^m(0)] /(2to) = 
Hf{0) and [.9ei(0) ~ 5sm(0)] /(2"i) = —Hg{0), the derivative boundary conditions 
appropriate for (5- functions, (iii) Incoming (incident), reflected and transmitted wave 
directions are defined by their group velocities, i.e., an electron incident from the left 
is transmitted with v > and reflected with v <0. 

If an electron incident on the interface from the SM with energy uj > A has wave 
vector q^, the four outgoing channels, with probabilities A, B, C, D, have respectively 
wave vectors q~ , —q'^, —k~, as shown in Fig. 11.11 C is the probability of trans- 
mission through the interface with a wave vector on the same (i.e., forward) side of 
the Fermi surface as q'^ (i.e., (7+ — > fc"*", not — fc~), whereas D gives the probability of 
transmission on the back side of the Fermi surface (i.e., q'^ — ?► —k~). B is the probabil- 
ity of intra-branch reflection, whereas A is the probability of Andreev (cross-branch) 
reflection. The steady state solution of system (I1.12p is 

I'Sm(z) = «'i„c(2) + «'rofi(2), ^'ei(2) = l'tians(z), 

where 

^,.Uz) = cI''Ac''^^^ + dhy-^''-^. (1.40) 
\voJ \uoJ 

Applying the boundary conditions, we obtain a system of four linear equations in the 
four unknowns a, b, c, and d, which wc solve at a fixed value for w. We introduce the 
dimcnsionless barrier strength Z — H/vp, and approximate = k~ = q~^ = q~ kp, 
on the basis that the ratio 2mA/A:p is small. The quantities A, B. C, D, are the ratios 
of the probability current densities of the specific transmission or reflection channels 
to the current of the incident particle, e.g. A = \JA/Jinc\, and so on. The conservation 
of probability requires that 

A + B + C + D = 1. (1.41) 

This result is useful in simplifying expressions for energies below the gap, oj < A, 
where there can be no transmitted electrons, so that C = D = 0. Then, Eq. p.4ip 
reduces simply to A — 1 — B. 



vo = 



1 - 



(a;2 _ A2)i/2 



(1.38) 



Andreev reflection at the interface between excitonic insulator and semimetal 23 



Table 1.2 Transmission and reflection coefficients for the SM-EI junction. A 
gives tlie probability of Andreev reflection (cross-branch), B of ordinary reflec- 
tion, C of transmission without branch crossing, and D of cross-branch trans- 
mission. Here 9 ^ cj^ + AZ'^cj^ + (1 + 'iZ*){A'^ - cj^) - 8Z^uj {A^ - uj'^Y^' , 
-f^Z'^ [vl - lig) + {iZ + 1/2) 2ul, and itg = 1 - = 1/2 [1 + {uj'^ - A'^y^'^/oj]. 

A B C D 

7^ 1 

No condensate i+z'^ ^ 

General form 

a;<A ^ 1- A 



tT FrP pyP FrF 



w > A 
No barrier 

cj < A 1 

a; > A 4 



Strong barrier 



"0 



^ ^ ^-Ki-t-j^)^ Z-Hu-„-v'f,) Z^uj-v-) 



We find 



7 



7 

uo (1 + i^) 
c = , 

7 

d = -^, (1.42) 

7 

7 = Z2 («2_y2) + (iZ+l/2)2u2. 

The probability coefficients are actually the currents, measured in units of v-p. For 
example, A = \ Ja\ /v-p = |a|^, and D ~ \df' / \vq — Uq\. The expression for the energy 
dependences of A, B, C, and D can be conveniently written in terms of the so-called 
coherence factors uq and vq. The results are given in Table [TT^ For convenience, in 
addition to the general results we also list the limiting forms of the results for zero 
barrier (Z = 0) and for a strong barrier [Z'^{uq — Vq) ^ 1], a.s well as for A = (the 
semimetal case). 

In the absence of disorder (Z ~ 0), the dependence of the transmission coefficient 
C{uj) on energies close to the gap is 



C{io) = 2^2 (w- A) /A w « A, (1.43) 

whereas below the gap the electron is totally Andreev-reflected and the transmission is 
zero. The ordinary reflection channel is completely suppressed {B = 0) as well as cross- 
branch transmission {D — 0). The situation is depicted in Fig. 11.5( a). Even above the 
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Energy co 

Fig. 1.5 Transmission and reflection coeiJicients at tlie SM-EI boundary vs quasiparticle 
energy uj. (a) Z — 0. (b) Z = —1. A gives tiie probability of Andreev reflection, B gives 
the probability of ordinary reflection, C gives the transmission probability without branch 
crossing, and D gives the probability of transmission with branch crossing. The parameter Z 
measures the interface transparency. 

gap, w > A, there is a high probabihty for Andreev reflection, which strongly depends 
on ui. For energies close to the gap, w w A, Andreev reflection is almost certain, Awl. 
This is the cause for the low value of interface conductance. The effect is washed out 
by the opacity of the interface: as \Z\ increases [Z — —1 in Fig. ll.Sf b)]. the total 
reflection probability A + B loses its dependence on w, and the dominant reflection 
channel changes from the Andreev one (A) into the ordinary one (B). 

The interface opacity Z is the handle to tune the effect of excitonic coherence 
on transport, as discussed in the next sections. Remarkably, the overall set of results 
of Table 11.21 is formally identical to the analogous quantities obtained for the N-S 
junction (e.g. compare with Table II of Ref. [22]), the only slight difference being the 
behavior for Z 7^ 0. In fact, due to the different boundary conditions, the coefficients 
of the N-S junction are even functions of Z, whereas those of Table FOl do not have a 
definite parity with respect to the sign of Z for w < A (there is a mistake in the entry 
corresponding to the sub-gap strong-barrier case appearing in Table II of Ref. [22]). 
Nevertheless, the expressions for coefficients in the strong-barrier case are the same 
for the SM-EI and the N-S junction. 

The appearance of coherence factors u and v in the coefficients of the SM-EI 
junction demonstrates that the electron-hole condensate strongly affects the transport 
and in general the wave function of carriers, by means of both inducing coherence on 
the SM side and altering transmission features. 

1.4.3 Proximity effect 

From the results of the previous two subsections it follows that the condensate on the 
right hand side of the junction induces pairing correlations in the normal phase on the 
left hand side, even if there the order parameter A is zero as interactions are switched 
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off. In tlie N-S junction tlie pairing induced by tlie proximity effect correlates electrons 
with opposite spins, whereas in the SM-EI junction it correlates electron-hole pairs. 

To show that the exciton condensate induces exciton order on the SM side it is 
sufficient to compute the built-in coherence (\l/o|V'b('") V'i('")l^o)- This would be zero 
in an isolated SM but by the proximity with the condensate acquires the value 



2 dujN{uj) cos 



k 

arctan | — | -|- 2 — z 



(1.44) 



with N{uj) being the density of states in the SM. Inside the gap (a; w 0) each quasi- 
particle contributes to the sum (|1.44p with a term ^ exp [i arctan (A/w) + 2\ijjz/vY\. 
The only coordinate dependence enters this expression via the phase factor, 2ijjz/v-p, 
which represents the relative phase shift of conduction- and valence-band components 
of the wave function. If w = 0, then these components keep constant relative phase 
arctan (A/cj) all the way to z = — oo, where no pairing interactions exist. Therefore, 
the reflected electron — equivalent to an incoming hole — has exactly the same veloc- 
ity as the incident particle, and will thus retrace exactly the same path all the way to 
z = — oo. At finite energy, the z dependent oscillations provide destructive interference 
on the pair coherence. Hence, the paths of incident and reflected electrons part ways 
away from the interface. Analogous considerations apply to the incident electron and 
to the Andreev-reflected hole in a sub-gap scattering event at the N-S interface [287]. 

1.5 A perfect insulator 

From the results for transmission and reflection probabilities obtained in Sec. 11.4.21 
derive in the linear response regime the values of the electrical and thermal conduc- 
tances of the SM-EI junction, G and Gt, respectively. The derivation is standard and 
it proceeds along the lines explained e.g. in Refs. P^53II215| . The Seebeck coefficient is 
zero due to the symmetry of the model |293il215j . Then, except for an additive phonon 
contribution to the thermal conductance, the interface thermoelectric properties are 
completely determined by G and Gt- 

Both G and Gt have an activation threshold at low temperature, T, proportional 
to the gap A, as shown in Fig. 11.61 (cf. the curves for the transparent barrier with 
Z = Q). At first sight, the transport properties of the SM-EI junction seem dramat- 
ically different from those of the N-S junction, as the latter may sustain an electric 
supercurrent at vanishing bias voltage whereas the former exhibits insulating behavior. 
However, a closer examination shows that the two junctions share essential features. 

In fact, the functional dependence of Gt on T and A is the same for both the 
SM-EI and N-S junctions. Remarkably, the time rate of entropy production, S, is 
the same very low value in both cases, pointing to the disipationless character of the 
flow of charge and heat. This is seen by the relation between S and the transport 
coefficients |293| . 

S = G{5Vf /T + AGT{5Tf /T'^ , (1.45) 

with 5T and 5V being respectively the temperature and voltage drops at the interface 
and A being the cross-sectional area. In the N-S junction the superfluid component 
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Fig. 1.6 (a) Electrical conductance G of the SM-EI junction vs temperature T. The curves 
correspond to different values of the dimensionless barrier opacity, Z = 0,1,3,5,7,9. For 
each value of Z, G is divided by the corresponding transmission coefficient when A = 0, i.e., 
Ca=o = {Z'^ + Kb is the Boltzmann constant, A is the interface cross-sectional area, 

N{eF) is the density of states evaluated at the Fermi energy ep. (b) Thermal conductance 
Gt of the SM-EI junction vs T. The curves correspond to Z = 0, 1, 3, 5, 7. 

does not carry any entropy. Therefore, the terms proportional to G and Gt only 
include the contribution of quasiparticles which, when they cross the N-S interface, 
experience the same electric and thermal resistance as electrons do across the SM-EI 
boundary. 

To shed light on the dissipationless motion of electrons in the linear transport 
regime, we vary the opacity of the SM-EI junction. The coherence between the two 
sides of the interface is diminished as the dimensionless barrier strength Z increases 
from zero (clean junction) to finite values (tunneling regime). Figure [T6l displays G 
and Gt for increasing values of Z. Since the transmission coefficient C{uj) decreases 
uniformly in the absence of any electron-hole pairing (cf. Table [TT^ . Ga=o = {Z^ + 
1)~^, we rescale conductances dividing them by Ca=o- Naively, we would expect that 
the insertion of an insulating layer would reduce the conductances. On the contrary, the 
effect is just the opposite: as Z increases, G/Ca=o a-nd G't/Ca=o increase, eventually 
reaching saturation in the tunneling regime. This shows that the exciton order induced 
in the SM side by EI makes the junction less conductive for charge and heat transport. 
The plot of the differential conductance {dl /dV)/CA=o vs the bias voltage V at low T 
fFig. 11.7]) allows clear monitoring of the transition from the transparent to the opaque 
limit, where transport is recovered. The effect is maximum for eV sa A and as T 0, 
when the differential conductance becomes proportional to C{eV) + D(eV). 
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Fig. 1.7 Differential conductance [dl / dV) / C a=o of tiie SM-El junction vs bias voltage V, 
computed at K^T j = 0.1. The curves correspond to different values of the dimensionless 
barrier opacity, Z = 0,1,3,5,7,9,11. For each value of Z, dl/dV is divided by the corre- 
sponding transmission coefficient when A = 0, i.e., Ca=o = {Z^ + 

1.5.1 Charge versus exciton current 

As anticipated at the end of Sec. 11.4.11 the transport features discussed above which 
distinguish the excitonic insulator from the normal insulating state may be explained 
by two alternate physical pictures. The conventional view is that electrons below the 
energy gap cannot contribute to transport as they are back scattered by the gap barrier, 
A, formed by the proximity effect of the EI. The less conventional view is to make use 
of the analogy with the N-S junction. Instead of counting the electrons in the valence 
band as negatively charged carriers of the current, we may start with the state with 
the valence band filled to the top as carrying zero current even under an electrical or 
thermal current and regard each unoccupied state in the valence band as a positively 
charged carrier — a hole — moving in the direction opposite to the electron. Then the 
reflected electrons are replaced by incoming holes towards the barrier. Therefore, the 
incident conduction electron and the valence hole may be viewed as a correlated pair 
moving towards the interface [Fig. I1.4f a')]. The novelty is that a constant electron- 
hole current moves from the SM to the EI below the gap, where electric transport 
is blocked. As the electron-hole pair approaches the interface from the SM side, the 
exciton current is converted into the condensate supercurrent: the global effect is that 
in the steady state an exciton current exists flowing constantly and reversibly all the 
way from the SM to the EI without any form of dissipation. 

The above scenario follows from the continuity equation for the electron-hole cur- 
rent. The probability density pe-h(T',i) for finding either a conduction- band electron 
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or a valence-band hole at a particular time and place is po-ii('",i) = |/| + 1 ^ |.9| • 
Thus, the associated continuity equation is 

+ V • Jc-h = 0, Jc-h = Jpair + «/cond- (1.46) 

One component of the electron- hole current, 

Jpair = ^Im{/*V/ + g*Wg} , (1.47) 

is the density current of the electron-hole pair, similar to the standard particle carrier 
J = m^^Im{/* V/— g*V.g} with an important difference in sign. The other component, 

V- Jeond = -4Im{r5A}, (1.48) 

depends explicitly on the built-in coherence of the electron-hole condensate A, and 
may be described as the exciton supercurrent of the EI state. 

Let us go back to our picture of A(z) smoothly varying in space (Fig. II. 2|) . with 
the junction being divided into small neighborhoods at position r and each being a 
homogeneous system. If < | Ao|, each electron wave function, solution of eqn (jl.28p . 
carries zero total electric current eJ, which is the sum of the equal and opposite 
incident and reflected fluxes, and finite and constant electron-hole current J^-h = 2vpn, 
with n a unit vector. When z — ^ — oo. far from the interface on the SM side, the 
supercurrent contribution Jcond is zero. As z increases and A(z) gradually rises, both 
J and Jo-h conserve their constant value, independent of z, since quasiparticle states 
are stationary. However, their analysis in terms of incident and reflected quasi-particles 
is qualitatively different. 

From the electron point of view, we see in Fig. ll.2r a) that the incoming conduction- 
band particle approaching the EI boundary sees its group velocity progressively re- 
duced (from time step 1 to time step 3), up to the classical turning point (time step 4) 
where it changes direction and branch of the spectrum: there is no net electric current. 

From the exciton point of view, as the contribution to the electron-hole current 
Jpair vanishes approaching the boundary, since the group velocity goes to zero at the 
classical turning point where the wave function becomes evanescent, Jpair is converted 
into the supercurrent Jcond. Excitons therefore can flow into the EI side without any 
resistance, and the sum Jc-h of the two contributions, Jpair and Jcond, is constant 
through all the space [Fig. I1.4f a)]. 

As cu exceeds |Ao|, J acquires a finite value and Je_h monotonously decreases. 
However, close to the gap, electron transmission to the EI side is still inhibited [cf. eqn 
(|1.43p ] by the pairing between electrons and holes of the condensate: an electron can 
stand alone and carry current only after its parent exciton has been "ionized" by 
injecting — say — a conduction-band electron or by filling a valence-band hole in the 
EI. The ionization costs an amount of energy of the order of the binding energy of 
the exciton, |Ao|. Therefore, as long as w ss |Ao|, the competition between exciton 
and electron flow favors Andreev reflection, which is the source of both electric and 
thermal resistances. 
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In equilibrium, there is no net charge or heat flow, since quasi-particles with v and 
—V compensate each other. However, if a heat current flows, the net drift velocity of 
electrons and holes locally "drags" the exciton supcrcurrent, which otherwise would 
be pinned by various scattering sources |93j . 

1.5.2 A concrete example 

As a concrete example of the aforementioned conversion of frcc-cxciton current into 
condensate supcrcurrent, consider the quasiparticlc steady state of Eq. (jl.40p . For 
< A, and in the EI (z > 0) have small imaginary parts which lead to an 
exponential decay on a length scale A, where 

The quasiparticles penetrate a depth A before the electron-hole current Jpair is con- 
verted to a supcrcurrent Jcond carried by the condensate; right at the gap edge the 
length diverges. For clarity, we define C and D here as the transmission probabilities 
at z 3> A, while for co > A plane-wave currents arc spatially uniform and we need not 
specify the position at which they are evaluated. 

When the interface is transparent, Z = 0, the steady state (|1.40p is specified by 
b = d = 0, a = vq/uq, and c ~ I/uq. Below the gap coherence factors uq and vq 
are complex and equal in modulus. For w < A, \a\ = 1, which means the incident 
conduction-band electron is totally reflected into the SM valence band. Thus, the 
electron- hole current Jpair carried in the semimetal equals 2vf, but Jpair of the cxcitonic 
insulator is exponentially small for z ^ 0. Explicitly, 



I c| 9 2 

Jpair = (IwoT + l^'o| ) Im 

m 



Letting fc+ ss fcp + i/(2A), we have 

Jpair = 2wFC-"/\ (1.50) 

The "disappearing" electron-hole current reappears as exciton current carried by the 
condensate. Recalling the definition of Jcond, 

9Jcond/5z = -4Im{rgA}, 



by integration we obtain 

Jeorrd = -4A |c|' ^ d z'e-^'/^ MKvo] = 2vF (l - e-^/^) . (1.51) 

This is the desired result, explicitly showing the supercurrent Jcond increasing to an 
asymptotic value as z — )■ oo, at the same rate as the free electron- hole current Jpair 
dies away. 
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1.6 Josephson oscillations between exciton condensates in 
electrostatic traps 

The Josephson effect is a a macroscopic coherent phenomenon which has been ob- 
served in systems as diverse as superconductors [T7|, superfluid Hehum p^ll93[[249] . 
and Bose-Einstein condensates of trapped ultra-cold atomic gases [411l234ll4lll54| . Since 
Josephson oscillations appear naturally when two spatially separated macroscopic wave 
functions are weakly coupled, they have been predicted for bosonic excitations in solids 
as well, like polaritons |2801I218| and excitons optically pumped in suitable semicon- 
ductor hetcrostructures [231) . However, unlike the polaritons, which have a photonic 
component allowing for easy detection, excitons stay dark unless they recombine ra- 
diatively. So far, it is unclear how the exciton Josephson effect could be observed. One 
proposal relied on the observation of quasiparticlc excitations from the spectral prop- 
erties of the emitted light |218j . A drawback of this idea is that spectral properties 
are not unambiguously linked to the condensed phase. In this Section we propose the 
interference of emitted photons as a direct probe of the exciton Josephson effect. 

Condensed excitons are predicted to act as coherent light sources |190| [68 lfT89] (see 
discussion below). If Josephson oscillations occur between two exciton traps, in princi- 
ple they can be probed by measuring the interference of the beams separately emitted 
from the traps. However, in the time interval before recombination and contrary to 
the polariton case, there are too few photons emitted for an adequate signal to noise 
ratio, so one has to average the signal over many replicas of the same experiment [103] . 
Here we show that such ensemble averaging blurs the signature of the Josephson effect 
except in the relevant case of exciton "plasma" oscillations |198| . For the latter the 
dipole energy difference between the traps modulates the visibility a of interference 
fringes, providing a means for detection. 

The subsections below are organized as follows: we first introduce the double quan- 
tum well system illustrating a feasible scheme to manipulate electrically the exciton 
phase (Sec. ll.6."T|) : we then set the theoretical framework (Sec. ll.6."2|) . and we eventually 
propose a correlated photon counting experiment (Sec. I1.6.3p . 

1.6.1 Electrical control of the exciton phase 

Consider a double quantum well where electrons and holes arc separately confined 
in the two layers. In experiments aiming at Bosc-Einstcin condensation of excitons, 
electron-hole pairs are optically generated at energies higher than the band gap, left 
to thermalize, form excitons, and, at sufficiently low temperature and high density, 
condense before radiative decay |112[|103] . Let z be the growth axis of the two wells 
separated by distance d. The electrons in the conduction band and holes in the valence 
band move in the planes z = d,0, respectively (Fig. ll.8( ). In the experiments [3T|l32| . an 
electric field is applied along z to suppress inter-layer tunneling, thereby quenching 
the exciton recombination. 

Fabrications [571I1061IH] of electrostatic traps with suitably located electrodes to 
provide lateral confinement for the excitons have been implemented. The double quan- 
tum well is sandwiched between two spacer layers, providing insulation from planar 
electrodes lithographed on both sides of the coupled structure. Each electrode controls 
a tunable gate voltage, Vg{x, y, z), which localizes in a region of the xy plane the field 
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Fig. 1.8 (a) Double quantum well energy profile along the growth direction z. Electrons and 
hole move freely in the xy planes at z = d and z = Q, respectively. Fz is the component of 
the electric field parallel to the exciton electric dipole. (b) Effective exciton potential profile 
vs X for the double-trap setup. 

component along z, Fz{x,y,z) ~ —dVg{x,y, z)/dz, while F^,- and Fy are small and 
can be neglected as well as the the dependence of on z. The vertical field Fz{x,y) 
makes the electrostatic potential energy of the exciton dipole depend on the lateral 
position, Ux{x,y) — —edFz{x,y) (e < 0) [cf. Fig. ll.Sr b)]. In this way, potential traps 
for excitons are designed with great flexibily, with in situ control of the height, width, 
and shape of the potential barriers [97 (11 04j . 

First, we focus on the quasi-equilibrium situation before radiative recombination, 
where excitons condense in two coupled electrostatic traps, both within the condensate 
coherence length. Figure fLST b) depicts the exciton potential profile Ux{x,y = 0) 
along the x axis, with a link between two identical traps. The potential barrier allows 
tunneling between the condensates in the two traps, whose macroscopic wave functions 
are 2i(a;, y, t) and S2(a;, y, t), respectively. The optical coherence in a single trap is of 
the form 

Eix, y, t) = y, 0, t)^b{x, y, d, t)) , (1.52) 

where (...) denotes quantum and thermal average in the grand canonical formalism. 
In eqn (|1.52p ^\[x,y,0,t) and 'ifb{x,y,d,t) are respectively the hole and electron 
destruction operators, with the vacuum being the semiconductor ground state with 
no excitons. With respect to definition (|1.13p . here by putting {x,y) = {x',y') we 
ignore the internal structure of the exciton relative-motion wave function. The reason 
is that we focus on the ideal BEC limit of dilute weakly interacting excitons, n a| ^ 1, 
with Ob being the two-dimensional effective Bohr radius and n the exciton density. 
Therefore, S(a;, y, t) is the macroscopic wave function for the ccnter-of-mass motion of 
excitons, which may be written in the form 
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(1.53) 



with Us being the density of the exciton condensate and ip the phase. We recall that 
only the relative phase between the two condensates has measurable effects [TTlllSOj . 
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For a gauge transformation of the gate potential Vg ^ Vg ~ c ^dx{t)/dt, which 
leaves the field F^^ unaltered, the field operators ^' gain a phase, 



*a ^ *a exp 
exp 



le , , 
he 

Tic 



(1.54) 



Throughout this section we indicate explicitly the reduced Planck's constant Ti. The 
macroscopic wave function, by eqn (jl.52p . also gains a phase, 

V^^+Yc [X{z ^d,t)- x{z = 0, t)] . (1.55) 

Hence, the frequency of time oscillation of the condensate is given by the electrostatic 
energy of the exciton dipole in the external field [T5], U = —edF^: 

ip = (^(0) + ledFJ , (1.56) 
n 

with if^^^ being the time-independent zero-field value. In the absence of the bilayer 
separation of the electrons and the holes, their gauge phases gained in the electric 
field would cancel each other resulting in no time dependence driven by U. Equation 
(|1.56p shows that the experimentally controllable dipole energy difference between the 
two traps depicted in Fig. ll.Sf b). AU = —ed{Fzi — -^22), drives the relative phase 
between the two condensates, thereby creating Josephson oscillations as a means for 
measuring the Josephson tunnel between the traps. 

1.6.2 Exciton Josephson oscillations 

We next introduce the usual two-mode description of inter-trap dynamics based on 
the Gross-Pitaevskn (GP) equation [2801 [211 [2311 [2Ml[2Q3[28i [Tig Exciton-exciton 
correlation [191j beyond the GP mean field may be neglected due to the repulsive 
character of the dipolar interaction between excitons in coupled quantum wells. The 
condensate total wave function solution is 

S(a:, y, t) = E,ix, y, iVi) e'"^^ + E^ix, y, N2) e*'^^ , (1.57) 

where both the trap population Ni (t) and the condensate phase ipi {t) possess the entire 
time dependence for the ith trap (i = 1,2), and Sj;(j;, y, Ni) is a real quantity, with 

dxJdyE^{x,y,N,)=N,it). (1.58) 

The dynamics of the GP macroscopic wave function E!(x, y, t) depends entirely on the 
temporal evolution of two variables, the population imbalance k{t) = {Ni — N2)/2 and 
the relative phase (j){t) ~ (pi — (p2 of the two condensates. Here we consider a time 
interval much shorter than the exciton lifetime (10 — 100 ns) and ignore the spin 
structure. Therefore, the total population is approximately constant, Ni{t) + N2{t) = 
N. 
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The equations of motion for the canonicahy coniugatcd variables hk and 4> are 
derived from the effective Hamiltonian 

Hj ^ E^— + AUk - y ViV2 - 4fc2 cos , (1.59) 

under the condition fc <C (Ref. [198] ). Ec — 2dfii/dNi is the exciton "charging" 
energy of one trap, where is the chemical potential of trap 1, whereas Sj is the 
Bardeen single-particle tunnelling energy, 

where m is the exciton mass. The single-particle orbital S,i{x,y) is defined through 
Sj(2;,y) = y/Nl^i{x,y). 

The various dynamical regimes associated to certain intitial conditions (fc(0), 0(0)), 
including tt oscillations and macroscopic quantum self-trapping, are exhaustively dis- 
cussed in Refs. [23611203] . Two cases are specially relevant: 

AC Josephson effect. Under the conditions AU » NEc/2, AU 6j, one easily 
obtains 

<j>{t) = -^t + m , fc=^sin</). (1.61) 

Equation (|1.61|) shows that, analogous to the case of two superconductors separated 
by a thin barrier, if the phase difference </> between the condensates is not a multiple of 
TT, an exciton supercurrent 2k flows across the barrier. Remarkably, in the presence of 
an electric field gradient along z, an exciton fiux oscillates back and forth between the 
two traps, with frequency AU /h. As an exciton goes through the barrier, it exchanges 
with the field the dipole energy acquired or lost in the tunneling process. The analogy 
with the AC Josephson effect for superconductors is clear: in that case a bias voltage 
V is applied across the junction, and the energy 2eV is exchanged between field and 
Cooper pairs, as the latter experience a potential difference of V when penetrating the 
potential barrier. 

Plasma oscillations. This case concerns small oscillations around the equilibrium po- 
sition {k, 0)cq ~ (0, 0). The Hamiltonian (|1.59p may then be linearized into the form 

Hj^^ (2^ + E,^ + hjN^' + Af/fc _ M . (1.62) 

It follows that both k and (j) oscillate in time with plasma frequency 



ujj^^^5j{NE,/2 + Sj). (1.63) 
Note that AU displaces the equilibrium position from (fc, (f))cq = (0, 0) to 



(1.60) 



(fc,</))eq = i-AUN6j/2{hojjf,0). 



(1.64) 
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detecto 




Fig. 1.9 Proposed experimental setup to measure the time correlation of the photons emitted 
from the two traps. The delay time r of one optical path is externally controlled. 

1.6.3 Correlated photon counting experiment 



Figure 11.91 illustrates the correlated photon counting setup which we propose to probe 
Josephson oscillations. The detector measures the intensity /(r) of the sum of the two 
beams separately emitted from the traps. A delay time t is induced in one of the two 
beams, as in Ref. |282j . The fields are simply proportional to the order parameters 
of the traps. In fact, S(j;, y, t) is associated with a macroscopic electric dipole moment, 



which couples to photons: 



P{t)=xP,{t)±iyPy{t), 



Px{t) ^ J dxdyx E{x,y,t), 



(1.65) 



(1.66) 



and similarly for Py. The built-in dipole {P{t)) ^ oscillates with frequency (/x + 
Ex)/fi, where Ex is the optical gap minus the exciton binding energy, and /i accounts 
for exciton-exciton interaction |1901 [68 lll89j . This macroscopic oscillating dipole is 
equivalent to a noiseless current, which radiates a coherent field |82j . 
Therefore, the measured intensity /(t) is 



/(r) = 2/o[l + (cos(/)(T))], 



(1.67) 



assuming that the fields emitted from the two traps have the same magnitude (and 
intensity /q) but different relative phase </>, which is evaluated at the delayed time r. 
/(r) may be written as 

/(r) =2/o[l + acos<?!)o(T)] , (1.68) 
where 4'q{t) is the phase averaged over many measurements, defined by the condition 

(sin[0(T)-<^o(T)]) =0, (1.69) 



and 
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a = (cos[0(r)-0o(r)]) (1.70) 
is the fringe visibility, i.e., the normalized peak-to- valley ratio of fringes, 

a = - , (1.71) 

-'max I -'min 

with /max (/min) being the maximum (minimum) value of /(t), and < a < 1. 

Equation (|1.68p has a few important caveats. Since /(t) is an average, the temporal 
inhomogeneous effect will blur the interference fringes, i.e., a < 1. Other dcphasing 
mechanisms include exciton recombination and inelastic exciton-phonon scattering 
[282], as well as inelastic |282j and elastic exciton-exciton scattering, which in first 
instance may all be neglected for short t, low T, and na| ^ 1, respectively. 

The most immediate caveat is that the exciton condensates in decoupled traps must 
acquire a relative phase if initially they condense separately without a definite phase 
relation. This scenario is analogous to the case of interference between independent 
laser sources first discussed by Glauber [82| and later studied experimentally for mat- 
ter waves [13] . Even though a one-shot measurement with sufficient resolution would 
display fringes, the relative phase 4>q{t) is also subject to intrinsic dephasing effects 
by quantum fiuctuations |82j . The latter are significant noise sources which affect a, 
when (j) and k are quantized into canonically conjugated quantum variables whereas 
in the GP theory used so far they were classical variables whose fluctuations where 
neglected. 

In the following, we quantize Hamiltonian (|1.59p in order to properly evaluate 
a = (cos [0(t) — 00 (■'■)]) as a quantum statistical average in finite traps. Therefore, we 
follow Ref. |197| and introduce the commutator 

<^,fcj = i. (1.72) 

The operator k now appearing in the quantized version of Hamiltonian (|1.59|) takes 
the form —id/d4>, whereas the ground state wave function is defined in the space of 
periodical functions of (f> with period 27r. If condensate oscillations are mainly coherent, 
the variance of (j> is small and the visibility is approximated by a = 1 — ^ ((A0)^). 

The most interesting case concerns plasma oscillations. For AU = 0, the ground 
state of the quantized version of the harmonic oscillator Hamiltonian (|1.62p is a Gaus- 
sian, with (po = 0, independent from r, and minimal spreading (A(/)^) « {Ec/2SjNy^^ . 
Therefore, the interferometer output is time-independent, / = 2/o(1 + Q!), showing con- 
structive interference, / > 2/o, with 



Not surprisingly, the visibility is controlled by the ratio Ec/SjN, reaching the maxi- 
mum a = 1 as Ec/6jN — )■ 0. In fact, a is given by the balance between the competing 
effects of tunnelling (oc 6jN), which enforces a well-defined inter-trap phase, and in- 
verse compressibility (oc Ec), which favors the formation of separate number states in 
the two traps, thus separating the two macroscopic wave functions. 
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1 2 3 4 5 

Delay time t ( 1 / Mj ) 



Fig. 1.10 (Color online). Beam intensity I{t)/Io vs delay time r, for AU/hcoj — 0.2, 0.5 
(dashed and solid lines, respectively), (a) T = and a = 1, 0.8 (black and red [light gray] 
lines, respectively), (b) a{T = 0) = 0.94 and IvbT/Hluj = 0, 1, 2 (black, red [light gray], and 
blue [dark gray] lines, respectively). 

A small finite value of AU in eqn (|1.62|) displaces the equilibrium position of the 
harmonic oscillator. Noticeably, the ground state is a coherent state with harmonic 
evolution of the average phase in time, 

0o(t) = sin (wjt), (1.74) 



whereas a is unchanged. The Gaussian probability density characteristic of the ground 
state for AU = now is simply carried back and forth in (j) space in the same motion 
as the expectation value 0o (''")• 

This key feature allows for directly monitoring r-dependent plasma oscillations of 
frequency wj through the photon correlation measurement (cf. Fig. 11.10]) . We evaluate 
the effect of thermal fluctuations on a via the formula 

E„ cxp {-{iEr\ 

where /3 = Xjk^T and 2(a„ - 1) = ((A^)^)^ is the variance of </> in the nth excited 
state whose energy is En- At low T, the excited states may be approximated as those 
of the harmonic oscillator, giving 



The above results are summarized by the formula 
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1 + a{T) cos I — ^ sin r 



(1.77) 



which is vahd for Ec/SjN <^ 1. 

For smah dipole energy variations, AU/Hluj <^ 1, the osciUating part within the 
square brackets of eqn (|1.77p may be written as —a{T) /2{AU/ hw j)^ siii^ ujj t. This 
shows that the visibihty a(T) of fringes, which oscihate hke sin^ wj r, is modulated 
by the experimentaUy tunable factor {AU/hujj)'^ /2. The dependence of /(r) on AU 
is illustrated in Fig. ll.lOl for two values of AU/fiUJj. As AU/fiUj is increased [from 
0.2 (dashed lines) to 0.5 (solid lines)], the amplitude of oscillations of I{t) shows a 
strong non-linear enhancement, providing a clear signature of Josephson oscillations. 
The oscillation amplitudes are larger for higher values of a [cf. Fig. ll.lOT a)]. and fairly 
robust against thermal smearing [cf. Fig. ll.lOT b)]. In fact. Fig. ll.lOT b) shows that the 
oscillation of /(r) is still clearly resolved for temperatures as high as T w hujj/kB- At 
even higher temperatures a{T) displays anharmonic effects |197| . with a{T) — as 
T-^00. 

The AC Josephson effect cannot be observed within our scheme. In fact, for large 
values of Af7, the term proportional to cos appearing in the Hamiltonian (|1.59|) may 
be neglected in first approximation, and the ground state wave function is a plane wave, 
(27r)^^/^ exp [m(/)], where fi is the integer closest to —AU/Ec. Since the probability 
density, (27r)~^, is constant, the phase is distributed randomly and the visibility is 
zero. Therefore, the correction to a coming from the inclusion in the calculation of the 
term neglected in (|1.59|) will be small and fragile against fluctuations. 



1.7 Conclusions 

The venerated topic of exciton Bose Einstein condensation is facing a rebirth as its 
investigation is fueled by advances in novel materials and technologies, as for double- 
layer semiconductors and graphene. Indeed, growing evidence shows that the concept 
of exciton condensation is a paradigm of many-body behavior. The observation of the 
coherence properties of the condensate, including superfluidity, is an important long- 
term goal of this field. In this Chapter we have explained how the exciton analogues of 
Andreev reflection and Josephson oscillations may be linked to measurable quantities. 
We hope that these ideas may stimulate further experiments along this path. 
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